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Abstract: The ratios of hydrogen Balmer emission line intensities in cataclysmic variables are signa- 
tures of the physical processes that produce them. To quantify those signatures relative to classifications 
of cataclysmic variable types, we applied the multivariate statistical analysis methods of prin cipal com- 
ponen ts analysis and discriminant function analysis to the spectroscopic emission data set o f lWilliamsl 
( 1983). The two analysis methods reveal two different sources of variation in the ratios of the emission 
lines. The source of variation seen in the principal components analysis was shown to be correlated 
with the binary orbital period. The source of variation seen in the discriminant function analysis was 
shown to be correlated with the equivalent width of the H/3 line. Comparison of the data scatterplot 
with scatterplots of theoretical models shows that Balmer line emission from T CrB systems is consis- 
tent with the photoionization of a surrounding nebula. Otherwise, models that we considered do not 
reproduce the wide range of Balmer decrements, including "inverted" decrements, seen in the data. 
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1 Introduction 

Cataclysmic variables (CVs) are semi-detached bina- 
ries in which a red (K or M) main-sequence, or slightly 
evolved, low-mass star loses mass to a white-dwarf pri- 
mary through overfilling its Roche lobe. They can be 
divided into two groups, the magnetic CVs (mCVs) 
and the non-magnetic CVs, according to whether the 
white-dwarf has a magnetic field that is strong enough 
to affect the accretion dynamics. Semi-detached bina- 
ries frequently possess an accretion disc, where gravi- 
tational energy is converted to kinetic energy and then 
to radiation. With mCVs the magnetic-field strength 
of the white-dwarf is high enough to disrupt the ac- 
cretion disc at the magnetosphere. The white-dwarf's 
magnetosphere will determine the inner radius of the 
accretion disc and the tidal force exerted by the sec- 
ondary star will truncate the outer region of the disc. 
If the magnetosphere is large enough, no disc will form. 

MCVs without an accretion disc are known as po- 
lars (or AM Herculis systems). In a polar the white- 
dwarf magnetic field is sufficiently strong (~ 10 MG) to 
also lock the white-dwarf and its companion star into 
synchronous (or almost synchronous) rotation with the 
binary orbit. The accretion flow near the white-dwarf 
is channelled by the field and a shock is often formed 
near the base of the the accretion column before the 
accreting gas settles onto the white-dwarf atmosphere. 
The shock heats the accreting gas, leading to the emis- 
sion of bremsstrahlung X-rays that can subsequently 
photoionize the cooler gas upstream. For systems con- 



sisting of a white-dwarf with a moderately strong mag- 
netic field (~ 1 MG), the white-dwarf spin and the 
orbital motion are not synchronous, although the sec- 
ondary may still be tidally locked and rotate with the 
orbit. These moderately magnetic systems are known 
as intermediate polars (IPs). Most IPs possess an ac- 
cretion disc, but the inner region of the disc is trun- 
cated by the magnetic field, and the flow is channelled 
in the white-dwarf magnetosphere, similar to polars. 
A shock can be formed in the accretion column, and 
the X-rays from the post-shock region can photoion- 
ize both the pre-shock gases in the accretion column 
and the cool gas in the accretion disc. It is worth 
noting that some IPs may not have a fully developed 
accretion disc. They probably have a thin accreting 
annulus, which couples loosely to the magnetosphere 
of the white-dwarf. 

Non-magnetic CVs form a heterogeneous class of 
objects, which include nova, dwarf nova and nova-like 
systems, primarily on the basis of their light curve be- 
haviour (although strict classification is sometimes not 
possible, e.g. some novae are also polars). They are 
asynchronously rotating systems, and they have a fully 
developed accretion disc with an inner edge extending 
to the white-dwarf surface. A substantial fraction of 
the accretion energy is dissipated in a transition layer 
between disc and the white-dwarf. The strong UV 
photons emitted from the transition region can be a 
main cause of photoionization in the accreting gas. 

CVs are strong emission line objects, particularly 
in quiescence. They all show a clear H Balmer emis- 
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sion sequences. He I, He II and high ionization lines of 
other species are also prominent in the optical spectra 
of many systems. The high ionization lines of the CVs 
are believed to be emission from photoionized gas in 
the accretion disc (for IPs and non- magnetic CVs), the 
accretion stream, or irradiatively heated surface of the 
companion star. Some of the lines may originate from 
hot spots in the system, for example in regions where 
the accretion stream impacts onto the accretion disc 
and where the accretion flow begins to couple with the 
magnetic field. In these cases the gas is often heated 
by shocks that are formed as consequence of an abrupt 
change in the dynamics of the accretion flow. (For re- 
views o f CVs and related accretion physics, see tffiarnerl 
11995T ) . ICropperl il99Cft . lwl feOOOft . lWu etafl J2003IU 

The lines from photoionized gases and shock heat- 
ed gases have different properties. An example is that 
of the relative strengths of the H Balmer emission 
lines. For in stance, simple pho toionization-recombina- 
tion models fOstcrbrock 1989) generally predict ratios 
of m 3.5 for Ha/H/3; yet observations show Ha/H/3 
which deviate significantly from these model predic- 
tions. Such deviation implies that the line formation 
process is more complex than simple photoionization 
and subsequent recombination. 

Here we analyse the H Balmer emission from CVs 
systematically using multivar iate statis t ical m ethods. 
We consider the systems in IWilliams! il983t) as the 
sample for our analysis. We search for correlations 
between the emission line ratios and CV classification 
and review the line formation mechanisms that may be 
responsible for the correlations. The paper is organized 
as follows. In JUwe review the multivariate statistical 
methods used to analyze the H Balmer emission data. 
These methods are discriminant function and princi- 
pal component analysis. In J|]we present the data in a 
four dimensional scatterplot space and give the results 
of the discriminant function and principal component 
analysis in that geometric space. We also review mod- 
els of Balmer emission relevant to CVs. In section 
we discuss possible physical causes, based on existing 
models, of the statistical behaviour of the data. Two 
significant directions of variation are found that corre- 
late with the binary orbital period and the equivalent 
width of the h/3 line. Our conclusion is given in ^ 



2 Statistical Methods 

We consider the systems listed in I Williamsl ljl983T) as 
our sample. The spectra were obtained in 1980 or 1981 
using the 1.3 m telescope at McGraw-Hill Observa- 
tory. As the data were collected with the same tele- 
scope/instruments and were reduced in a similar man- 
ner, this ensures that any variations present are less 
likely to be caused by different experimental settings 
and data analysis procedures. The emission line ratios 
Ha/H/3, H7/H/3, rlS/R/3 and He/H/3 were derived from 
the spectroscopic data presented in IWilliamsl <ll983T) 
and represent integrated line intensities above the con- 
tinuum. Only systems for which all four H Balmer 
emission lines were measured were retained and sources 



since i dentified as non-CVs or were not listed in lDownes et alJ 
(12001 (TU Leo, Z And, HK Sco and CL Sco) were 
discarded. The remaining data (representing 95 of 153 
original spectra) are reproduced in Table Q along with 
four classification schemes. None of the systems con- 
sidered here where in outburst when the data were 
taken. There are three observations of two AM CVn 
stars which have very little hydrogen. However, the 
AM CVn data were retained because their H Balmer 
line ratios were reported. We have also not given mul- 
tiple observations of a single star any special status in 
our analysis assuming, as a first approximation, that 
the data represent a random sample from the observ- 
able CVs. 

Of the classification schemes considered, the first 
classificatio n is that from the A tlas of CVs: The Liv- 
ing Edition ijDownes et al.ll200ll) as of February, 2005. 
Based on the Downes et al. classification, three other, 
simpler, classification schemes were produced. For clas- 
sification scheme Class 1, the sources were grouped 
into one of four groups: (i) dwarf novae, (ii) polar & 
IP (i.e. mCVs), (iii) nova & nova-like, and (iv) dou- 
ble degenerate. (Double degenerate systems contain 
two white-dwarfs instead of a white-dwarf and a red 
dwarf.) All sources classified as DQ Herculis systems 
were put in the polar & IP (mCV) group unless the 
DQ Her classification was listed as very uncertain. The 
results of a discriminant function analysis (see below) 
of the groups in Class 1 lead to the definition of the 
groups in Class 2. That is, the groups of Class 2 were 
determined by a statistical analysis of the physically 
distinct groups of Class 1. For classification scheme 
Class 2 the sources were split into two groups, (i) dwarf 
nova and (ii) other, based on the grouping of Class 1. 
For classification scheme Class 3 the sources were split 
into two groups, (i) dwarf nova and (ii) mCVs, based 
on the grouping of Class 1, which excluded the nova 
& nova- like and the double degenerate groups. The 
grouping of Class 3 was defined in an attempt to find 
statistical differences between objects that are similar 
but for the magnetic field of the primary. 

Discriminant functions were constructed, using the 
SPSS statistical software (SPSS Inc., Chicago, USA), 
to maximally separate the two groups of Class 2 and 
Class 3. The discriminant functions were of the form 

4 

Z = C + ^ CiXi (1) 

1=1 

where Xi are the emission line ratios and d are coeffi- 
cients such that the values of z for the centroids (mean) 
of the two groups are maximally different (co is arbi- 
trary). Specifically we used x\ to represent Ha/H/3, X2 
to represent H7/H/3, £3 to represent H5/H/3 and X4 to 
represent He/H/3. Discriminant functions may be used 
to quantify the separation of more than two groups; 
we originally computed a discriminant function for the 
classification scheme of Class 1 but found significant 
separations only for the groups subsequently defined 
in Class 2 and 3. The statistical significance of the dis- 
criminant function was characterized with Wilk's A, a 
statistic that generalizes Fisher's analysis of variance 
F statistic to the multivariate situation, and an associ- 
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Table 1 : Balmer line emission data as taken from Wil liams! l)l983fl together with the classi- 
fications used in the discriminant function analysis. 



Star" 


Ha 

H/3 


H7 
573 


H(S 

H/3 


He 

H/3 


ACV 6 
Class 


1 


Class c,d ' 
2 


3 


EW(/3) J 

(A) 


DW(P) y 
(km/s) 


^orb 
(days) 


V368 Aql 


0.98 


0.92 


0.87 


0.85 


na 


N&NL 


O 








0.3452 


V603 Aql 


0.49 


1.33 


1.50 


1.70 


na/dq:: 


N&NL 


O 




7.7 




0.1385 


V603 Aql 


0.44 


1.22 


1.43 


1.55 


na/dq:: 


N&NL 


O 




6.2 




0.1385 


DN Gem 


0.53 


1.17 


1.43 


1.58 


na 


N&NL 


O 




8.4 


562 


0.12785 


DQ Her 


0.79 


1.13 


1.35 


1.30 


na/dq 


P&IP 


O 


M 


33.3 


505 


0.193621 


DQ Her 


0.59 


1.09 


1.26 


1.30 


na/dq 


P&IP 


O 


M 


16.7 


519 


0.193621 


DQ Her 


0.59 


1.13 


1.21 


1.26 


na/dq 


P&IP 


O 


M 


20.5 


460 


0.193621 


V533 Her 


0.43 


1.29 


1.62 


1.64 


na 


N&NL 


O 




4.9 


396 


0.147 


DI Lac 


0.58 


1.13 


1.30 


1.22 


11a 


N&NL 


O 








0.543773 


HR Lyr 


0.52 


1.25 


1.37 


1.43 


11a 


N&NL 


O 










BT Mon 


0.75 


1.04 


1.08 


1.08 


11a 


N&NL 


O 




23.3 


343 


0.338814 


BT Mon 


0.92 


0.99 


0.96 


0.93 


11a 


N&NL 


O 




35.0 


592 


0.338814 


BT Mon 


0.87 


1.04 


1.08 


1.09 


11a 


N&NL 


O 




30.6 


644 


0.338814 


V841 Oph 


0.78 


1.05 


1.17 


1.13 


11b 


N&NL 


O 




1.6 




0.6014 


GK Per 


1.25 


0.71 


0.65 


0.62 


na/dq 


P&IP 


O 


M 


10.8 


309 


1.9968 


X Ser 


0.56 


1.06 


1.37 


1.17 


nb: 


N&NL 


O 




7.1 




1.48 


T CrB 


2.97 


0.65 


0.46 


0.29 


nra 


N&NL 


O 








227.5687 


RS Oph 


3.37 


0.60 


0.43 


0.43 


nra 


N&NL 


O 




12.7 




455.72 


TPyx 


0.56 


1.24 


1.41 


1.50 


nrb 


N&NL 


O 




9.7 


554 


0.076223 


V1017 Sgr 


1.42 


0.47 


0.50 


0.34 


nb 


N&NL 


O 








5.714 


RX And 


0.84 


1.16 


1.21 


1.30 


ugz 


DN 


DN 


DN 


58.3 


664 


0.209893 


RX And 


0.54 


1.52 


1.77 


1.79 


ugz 


DN 


DN 


DN 






0.209893 


AR And 


0.97 


1.26 


1.29 


1.63 


ug 


DN 


DN 


DN 


38.5 


702 


0.163 


AR And 


0.84 


1.16 


1.63 


1.71 


ug 


DN 


DN 


DN 


39.2 


379 


0.163 


UU Aql 


0.97 


1.05 


1.20 


1.24 


ug 


DN 


DN 


DN 


50.4 


454 


0.163532 


UU Aql 


0.97 


1.16 


1.39 


1.55 


ug 


DN 


DN 


DN 


80.9 


670 


0.163532 


SS Aur 


1.06 


1.20 


1.27 


1.50 


ug 


DN 


DN 


DN 






0.1828 


SS Aur 


1.13 


1.08 


1.16 


1.37 


ug 


DN 


DN 


DN 


122.0 


569 


0.1828 


FS Aur 


0.83 


1.12 


1.25 


1.45 


ug 


DN 


DN 


DN 


60.8 


903 


0.0595 


FS Aur 


0.86 


1.06 


1.14 


1.38 


ug 


DN 


DN 


DN 


35.3 


571 


0.0595 


Z Cam 


0.48 


1.17 


1.36 


1.56 


ugz 


DN 


DN 


DN 






0.289841 


HT Cas 


0.75 


1.16 


1.26 


1.64 


ugsu 


DN 


DN 


DN 


97.7 


943 


0.073647 


HT Cas 


0.78 


1.13 


1.38 


1.72 


ugsu 


DN 


DN 


DN 


117.0 


1032 


0.073647 


WW Cet 


0.96 


1.22 


1.46 


1.70 


ugz: 


DN 


DN 


DN 


42.7 


849 


0.1758 


WW Cet 


0.64 


1.13 


1.22 


1.27 


ugz: 


DN 


DN 


DN 


25.4 


673 


0.1758 


1E0643.0 


0.85 


1.08 


1.33 


1.59 


ug/ugz 


DN 


DN 


DN 


96.5 


776 


0.216778 


1E0643.0 


1.18 


0.91 


0.82 


1.14 


ug/ugz 


DN 


DN 


DN 


81.3 


584 


0.216778 


1E0643.0 


1.02 


0.72 


0.78 


0.86 


ug/ugz 


DN 


DN 


DN 


43.5 


161 


0.216778 


1E0643.0 


0.99 


0.96 


1.02 


1.22 


ug/ugz 


DN 


DN 


DN 


64.0 


506 


0.216778 


1E0643.0 


0.53 


1.25 


1.47 


1.67 


ug/ugz 


DN 


DN 


DN 


6.1 




0.216778 


1E0643.0 


0.81 


1.08 


1.13 


1.33 


ug/ugz 


DN 


DN 


DN 


49.1 


873 


0.216778 


SY Cnc 


0.73 


0.99 


1.09 


1.09 


ugz 


DN 


DN 


DN 


13.7 


498 


0.38 


SY Cnc 


0.54 


1.42 


1.56 


1.64 


ugz 


DN 


DN 


DN 






0.38 


YZ Cnc 


0.72 


1.34 


1.45 


1.70 


ugsu 


DN 


DN 


DN 


80.1 


941 


0.0868 


YZ Cnc 


0.66 


1.33 


1.58 


2.06 


ugsu 


DN 


DN 


DN 


131.0 


906 


0.0868 


AK Cnc 


0.72 


1.00 


0.91 


1.01 


ugsu 


DN 


DN 


DN 






0.0651 


SS Cyg 


1.13 


0.89 


0.92 


1.05 


ugss 


DN 


DN 


DN 


69.5 


477 


0.27513 


SS Cyg 


1.05 


0.91 


0.86 


0.97 


ugss 


DN 


DN 


DN 


67.0 


513 


0.27513 


SS Cyg 


0.70 


1.18 


1.29 


1.46 


ugss 


DN 


DN 


DN 


30.2 


786 


0.27513 


EM Cyg 


0.60 


1.10 


1.33 


1.42 


ugz 


DN 


DN 


DN 


2.6 




0.290909 


AB Dra 


0.59 


1.16 


1.35 


1.50 


ugz 


DN 


DN 


DN 


19.6 


973 


0.152 


IR Gem 


0.77 


1.16 


1.37 


1.60 


ugsu 


DN 


DN 


DN 


80.6 


554 


0.0684 


AH Her 


0.97 


0.96 


0.98 


1.06 


ugz 


DN 


DN 


DN 


26.5 


682 


0.258116 


EX Hya 


0.49 


1.16 


1.34 


1.59 


ug/dq 


P&IP 


O 


M 


54.4 


939 


0.068234 


EX Hya 


0.55 


1.12 


1.25 


1.38 


ug/dq 


P&IP 


O 


M 


75.2 


1492 


0.68234 


EX Hya 


0.61 


1.18 


1.30 


1.39 


ug/dq 


P&IP 


O 


M 


80.7 


1144 


0.68234 
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Table 1: — Continued. 



Star" 


Ha 


H7 


US 


He 


ACV b 




Class cd 


e 




DW(/3) 9 


P u fl 
orb 




H/3 




H/3 


ii/j 


Class 


l 


2 


3 


(A) 


(km / s) 




T Leo 


0.60 


1.42 


1.50 


1.54 


ugsu 


DN 


DN 


DN 


115.0 


457 


U.UJOOZ 


X Leo 


1.05 


1 21 


1.38 


1.67 


ug 


DN 


DN 


DN 


22.3 


919 


0.1644 


CN Ori 


0.49 


1 33 


1.54 


1.59 


ugz 


DN 


DN 


DN 






n 1 R3i qq 


CZ Ori 


1.05 


1.10 


1.13 


1.37 


ug 


DN 


DN 


DN 


43.5 


503 


0.2189 




n fin 


1 21 


2 81 


O.OJ 


ugss 


DN 


DN 


DN 


A 7 




n 37 A(\ 


W7, So-p 
vv Zj oge 


0.47 


1 30 


1.55 


1.71 


ii(Twr7 tr\ri 
UtiVVZj/ t-H-J 


P&IP 


Q 


M 


27.8 


725 


U. UtJUUOO 


SW TTMa 

O VV vJ Ivlel 


0.42 


1.35 


1.51 


1.58 


1 1 <TC 11 / f \ C\ 

ugsu/aq 


P&IP 


o 


M 






U.UJUOIJ 


SW TIMa 


0.66 


1.22 


1.41 


1.63 


1 1 <TC 11 / A"1 C\ 

ugsu/uq 


P&IP 


o 


M 


64.3 


657 


u. uouo ± o 


SW TIMa 


0.68 


1.18 


1.39 


1.55 


1 1 ffc ii /nn 

ugsu/aq 


P&IP 


Q 


M 


71.6 


825 


U.UUUO J. (J 


TW Vir 


0.79 


1 21 


1.41 


1.50 


ug 


DN 


DN 


DN 


69.8 


805 




TW Vir 


0.91 


1.26 


1.38 


1.99 


u § 


DN 


DN 


DN 


96.5 


601 


189fS7 


TW Vir 


0.91 


1.21 


1.35 


1.75 


ug 


DN 


DN 


DN 


73.3 


680 


n i &9^7 


AE Aqr 


1.30 


71 


0.61 


0.52 


(In 

aq 


P&IP 


Q 


M 


21.7 


369 


n 4i 1 R^ifi 

U. 4: J. -LUtJU 


TT Ari 


0.44 


1 30 


1.47 


1.57 


yy/dq: 


P&IP 


Q 


M 


7.0 




1 3755 


G61— 29 


0.64 


1.21 


1.42 


1.50 




DD 


o 








(1 D3933Q 


AM rVn 


0.37 


1.43 


1 71 


1.83 


ibwd 


DD 


o 








fl 01 1 QD7 
u . u ± iyu i 


AM rVn 


0.37 


1.43 


1.75 


1.93 


ibwd 


DD 


Q 








01 1 QD7 
u . u ± iyu i 


YY Dra 


1 41 


1 OS 

1 .UO 


i ie 

_L. lo 


1 9Q 


Art 

aq 


P^TP 


n 


M 




73Q 

I O.J 


fl 1 R537A 


9 Arm 1 —997 


1.09 


0.96 


0.91 


0.97 




P&IP 


o 


M 


43.9 


901 




AM Her 

JT.1V1 llt.1 


0.98 


1.43 


1.59 


1.84 




P&IP 


Q 


M 


18.4 


230 


1 9KQ97 


AM FTpr 


0.85 


1 27 


1.39 


1.55 




P&IP 


Q 


M 


15.2 


232 


fl 1 9KQ97 


AM FTer 

JT.1V1 llt:l 


0.72 


1.26 


1.46 


1.59 




P&IP 


o 


M 


55.4 


517 


D 1 9KQ97 


AM FTer 

JT.1V1 lltJl 


1.01 


1.16 


1.34 


1.47 




P&IP 


o 


M 


48.8 


355 


fl 1 9KQ97 


V42fi Onh 

V Tt-^U Mil 


0.66 


1.18 


1.34 


1.45 




P&IP 


o 


M 


14.9 


813 




V442 Oph 


0.54 


1.10 


1.22 


1.25 


vy 


N&NL 


o 




4.8 


469 


0.12433 


V2051 Oph 


0.60 


0.98 


1.18 


1.35 


ugsu 


DN 


DN 


DN 


124.0 


853 


0.062428 


VV Pup 


0.62 


1.25 


1.24 


1.20 


am 


P&IP 


O 


M 


93.9 


513 




VSge 


0.50 


1.25 


1.59 


1.71 


nlv 


N&NL 







52.9 


739 


0.514197 


VSge 


0.61 


1.24 


1.43 


1.62 


nlv 


N&NL 







62.2 


645 


0.514197 


N Sgr 1962 


1.96 


0.85 


0.74 


0.72 


nra 


N&NL 







59.6 


208 




RW Tri 


0.59 


1.17 


1.26 


1.34 


UX 


N&NL 


o 








0.231883 


UX UMa 


0.47 


1.31 


1.51 


1.79 


UX 


N&NL 







1.1 


360 


0.196671 


AN UMa 


0.52 


1.14 


1.22 


1.24 


am 


P&IP 





M 


28.1 


783 


0.79753 


AN UMa 


0.51 


1.14 


1.19 


1.36 


am 


P&IP 


o 


M 


24.2 


661 


0.79753 


2H2215-086 


0.58 


1.21 


1.42 


1.42 


dq 


P&IP 





M 


25.0 


382 


0.20206 


2H2215-086 


0.47 
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a G61-29 = GP Com, YY Dra = DO Dra, 2A0311-227 = EF Eri, N Sgr 1962 = V3890 Sgr, 2H2215-086 
= FO Aqr, 2H2252-035 = AO Psc, 0633+ 71 = BZ Cam 

'Atlas of Cataclysmic Variables (ACV) IDownes et al.l200lT nomenclature: ug = U Gem variable (dwarf 
nova), ugz = U Gem variable (Z Cam subtype), ugss = U Gem variable (SS Cyg subtype), ugsu = U Gem 
variable (SU UMa subtype), ugwz = U Gem variable (WZ Sge subtype), na = Fast nova, nb = Slow nova, 
nra = Recurrent nova - giant donor, nrb = Recurrent nova - non-giant donor, nlv = Novalike variable (V 
Sge subtype), ux = Novalike variable (UX UMa subtype), vy = Novalike variable (VY Scl subtype - systems 
which undergo low states), am = AM Herculis variable (synchronous rotators), dq = DQ Herculis variable 
(non-synchronous rotators), ibwd = Interacting binary white-dwarf, : = Uncertain, :: = Very uncertain 

C DN = dwarf nova, P&IP = polar and IP, N&NL = nova and NL, DD = double degenerate 

^DN = dwarf nova, O = other 

e DN = dwarf nova, M = magnetic 

^EW = equivalent width 

9 DW = Doppler width 

^Orbital period data from Dowries et al. 12001) 
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Figure 1: Scatterplot of the emission line data given in Tabled IP = intermediate polar, NL = nova- like, 
DD = double degenerate (which are systems containing two white-dwarfs) . Individual systems are difficult 
to identify in general at this scale but four outlying nova and nova-likes and one outlying dwarf nova are 
obvious. From top to bottom in the H7/H/3 versus Ha/H/3 box, the nova and/or nova-likes are RS Oph, 
T CrB, N Sgr 1962 and V1017 Sgr. The outlying dwarf nova is RU Peg. See text for a discussion of those 
stars. 
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ated p- value that gives the probability of incorrectly re- 
jecting the null hypothesis of no difference between the 
two groups. The null hypothesis is generally rejected 
and statistically significant differences are said to be 
found when p < 0.05. The discriminant function test 
of two group separation on the basis of the multivariate 
line ratio data is mathematically equivalent, under the 
assumption of normal homoscedasticity (which means 
that all groups are assumed to be multivariate nor- 
mally distributed with the same covariance matrices), 
to Hotelling's T 2 test. The T 2 statistic tests if the 
mean vectors of two groups are significantly different 
and is equivalent to A because T 2 and A give th e same 
likelihood ratio test with two groups 
Histograms, along with the best fit normal distribu- 
tion, of the discriminant function values were plotted 
to show the group separation provided by the discrim- 
inant function. To aid in the interpretation of the re- 
sulting discriminant functions a principal component 
analysis (PCA) was also conducted. 

PCA finds the best description of the data in scat- 
terplot space assuming the y are a sample from a multi- 
variate normal population (Renchcr 2002). Surfaces of 
constant a (standard deviation from the mean vector) 
of a multivariate normal distribution form ellipsoids in 
scatterplot space. So PCA finds an ellipsoid in scat- 
terplot space such that the major axis of the ellipsoid 
lies in the direction of maximum variation of the data; 
this is the first principal component direction. Subse- 
quent principal component directions are directions of 
maximal variation that are geometrically orthogonal 
to the previous principal component directions. The 
ellipsoid is the geometrical expression of the data's co- 
variance matrix - the eigenvalues and eigenvectors of 
the covariance matrix define the size and direction of 
the ellipsoid. It should be noted that PCA consid- 
ers all the data as coming from a single multivariate 
normally distributed population with one mean vector 
while discriminant functions assume that each group 
comes from it's own population with different mean 
vectors (but the covariance matrices of each popula- 
tion are assumed equal). The discriminant function's 
level surfaces will generally contain the principal axes 
of the individual group ellipsoids (since they are the 
same, except translated, under the assumption of ho- 
moscedasticity) to achieve maximal group separation. 
From this it can be seen that the normal to the dis- 
criminant function's level surfaces should be roughly 
perpendicular to the PCA's first component direction 
defined from the ellipsoid that considers all the data 
together. 

For comparison to the data, scatterplots of H Bal- 
mer emission line ratios from accretion discs predicted 
by theoretical models were constructed. Four mod- 
els were consi dered. The first is the optically thin 
nebula model <|Osterbrocklll989ft . which assumes that 
the line formation is dominated by photoionization- 
recombination. The second model considered the tem- 
perature and density distribution i n the accretion disc 
expli citly via an a-disc prescription iShakura fc Sunvaevl 
1973). The line calculations were carried out assum- 
ing a local thermodynamic equilibrium (LTE) condi- 
tion with emission primarily from the outer region of 



(a) Dwarf novae 




Class 2 Discriminant Function Value 



(b) Other CVs 




Class 2 Discriminant Function Value 



Figure 2: Histograms of discriminant function val- 
ues for discriminating dwarf novae from all other 
CVs. All other CVs include nova in low states and 
mCVs. (a) Dwarf novae, (b) Other CVs including 
polars, IPs, nova remnants and nova-like variables. 
Overlain on each histogram is a normal curve hav- 
ing the same mean and standard deviation as the 
histogram. N = number of systems. 
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Figure 3: Histograms of discriminant function val- 
ues for discriminating dwarf novae from mCVs. (a) 
Dwarf novae, (b) MCVs including polars and IPs. 
Overlain on each histogram is a normal curve hav- 
ing the same mean and standard deviation as the 
histogram. N = number of systems. 



the d isc, which is opt ically thin to the continuum emis- 
sion <IWilliamsl ll980h The third model assumed that 
the emission region is isothermal and employed a non - 
LTE radiative transfer formulation (Wil liamsl I1991T) 
in the line calculations. Finally we co nsidered the 
chromospheric models of. Williamsl l|l995f) where a hot, 
photoionized chromosphere above the accretion disc is 
modelled. 

3 Analysis 

3.1 Observed properties of H lines 
in the CV samples 

The scatterplot of the emission line ratio data from 
Table Q is given in Fig. Q It may be discerned that 
the data tend to fall on a line in the 4D scatterplot 
space. Visible in the scatterplot are some obviously 
outlying systems. Those systems are as follows. RS 
Oph is a recurrent nova with past outbursts observed 
in 1898, 1933, 1958, 1967, 1985 and 2006, T CrB is a 
recurrent nova with past outbursts observed in 1866 
and 1946, and N Sgr 1962 (V3890 Sgr) is a recurrent 
nova with past outbursts observed in 1962 and 1990. 
RS Oph, T CrB and V3890 Sgr are all of the T CrB 
subclass of recurrent nova which have M giant secon- 
daries, P or b ^ lOOd and a high rate of mass accre- 
tion, Mi ^ lx 10 -8 Mq yr _1 , onto a white-dwarf pri- 
mary whose m ass is close to the Chandrasekhar mass 
llWarneii rr995h The re are no other T CrB type sys- 
tems in the data set. iDuerbeckl (Il987ft lists V1017 Sgr 
as a possible symbiotic star or recurrent nova. RU Peg 
is a U Gem type of dwarf nova with no other known 
characteristics that would suggest that it is different 
from other dwarf novae. 

The discriminant function that maximally split the 
Class 2 groups of dwarf novae versus others was 

z = 1. 799-0.152x1- 1.927x 2 -9.714x 3 +9.226x 4 . (2) 

The value of Wilk's lambda for this function was A = 
0.729 with p = 9x 10~ 6 . The unit vector V 2 = [ -0.011 
-0.142 -0.718 0.681 ] T is normal to the level surface 
hyperplanes of the discriminant function for Class 2. 
Histograms of the distribution of the values for the 
Class 2 discriminant function are given in Fig. [5] 

The discriminant function that maximally split the 
Class 3 groups of dwarf novae versus mCVs was 

z = -0.583 + 1. 125a;i - 1.229x2 - 7.812x 3 + 7.757x 4 . 

(3) 

The value of Wilk's lambda for this function was A = 
0.791 with p = 0.004. The unit vector V 3 = [ 0.101 
-0.110 -0.702 0.697 ] T is normal to the level surface 
hyperplanes of the discriminant function for Class 3. 
Histograms of the distribution of the values for the 
Class 3 discriminant function are given in Fig. |3 The 
angle between V2 and V3 is 6.8°. 

The direction first principal component (direction 
of maximum data variance) is given by the unit vector 
Vi = [ -0.441 0.510 0.530 0.514 ] T . The direction of a 
least squares line fit through the data in the 4D scat- 
terplot space is given by the unit vector V4 = [ —0.816 
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0.236 0.339 0.405 ] T and the angle between Vi and V 4 
is 29.8°. The angles between the first principal compo- 
nent direction, Vi, and the two discriminant function 
directions, V2 and V3, are 95.6° and 96.6° respec- 
tively. 



3.2 Line characteristics of theoreti- 
cal models 

3.2.1 Optically thin photoionized plasmas 

It is useful to consider emission from photoionized gas, 
which may be present in a disc, accretion stream or in 
surrounding nebulousity left over from previous nova 
eruptions since it is exp ected that all CVs experience 
periodic nova eruptions (IWarnerl H995'). For this pur- 



Dose the photoionization-recombination models in lOsterbrO' 
( 1989) are relevant and the resulting scatterplot points 
are given in Fig.Q] After excitation via boundary layer 
X-ray and UV illumination, optical recombination line 
emission may be modelled under two situations known 
as Case A, in which all emitted photons escape from 
the gas, and Case B, where the optical depth is high 
enough to convert every Lyman-line photon into lower- 
series photons plus either a Lyman a photon or two 
continuum photons through scattering. For Case B 
the assumed electron number density, JV C is relevant, 
and the models illustrated in Fig. 2] assume N c = 10 2 , 
10 4 and 10 6 cm -3 . It may be seen that the depen- 
dence of the photoionization caused Balmer line ratios 
on temperature and gas density is relatively weak but 
that significantly non-flat Balmer decrements are pre- 
dicted. In addition, the trend of scatterplot positions 
for higher density gas versus lower density gas for the 
photoionization models is in the same general direc- 
tion as the two disc models described below with, for 
example, the higher density cases being positioned in 
the upper left hand part of the scatterplot and the 
lower density discs in the lower right hand part for 
the Ha/H/3 versus H7/H/3 (abscissa versus ordinate) 
scatterplot. 

3.2.2 Accretion discs with optically thin 
continuum 



Williams ( 1980.) presented calculations of H Balmer 
lines from accretion discs, with an outer region opti- 
cally thick to the lines but optically thin to the con- 
tinuum. The emission lines arise directly from a hot 
disc instead of from a disc chromospheric surface layer 
heated by X-ray and UV radiation from the accretion 
boundary layer near the white-dwarf surface. An a- 
disc prescription is used to determine the radial den- 
sity and temperature profile, and the line-centre in- 
tensities are given b y the local Planc k function at the 
relevant frequencies. IWilliamsl <ll980t) calculations pre- 
dict that the position of a line in the emission line ratio 
scatterplot space depends on the disc mass flow rate 
M (Fig. |1J. The predicted equivalent width of the 
H/3 line, EW(/3), varies from 0.2 A for M = 10" 9 
M Q yr" 1 to 65 A for M = 10~ 12 M Q yr" 1 . The 
cases represented in Fig. 0] are for a disc of outer ra- 
dius r<j = 4 x 10 10 cm and a white-dwarf of radius 



Ri — 6 x 10 s cm and mass M\ = 1.0 Mq . The assumed 
disc radius r d corresponds to 1/3 the separation of two 
1 Mq stars having an orbital period of 4 hr. 

We note that in order to have the assumption of an 
optically thin continuum satisfied, the accretion rate 
of the system must not greatly exceed 10 -10 Moyr -1 . 
Most CVs are, however, observed to have higher accre- 
tion rates, so that models assuming an optically thin 
continuum are not always applicable. 

3.2.3 Non-LTE emission from isothermal 
medium 

The pred i ctions of a non-LTE model as computed by 
IWilliarnf] l)l99lfl are summarized in Fig. The model 
ssumes an isothermal horizontally infinite disc llWilliams fc Shipmanl 

with a Gaussian dropoff in density with vertical 
height from the midplane. The H number density No 
(in cm -3 ) represents the value in the disc midplane, 
which is generally a monotonic function of disc mass 
flow rate. The position of the line ratios on the scat- 
terplot are more influenced by the value of No than by 
the temperature with the higher density discs being 
represented in the upper left hand part of the scatter- 
plot and the lower density discs being represented in 
the lower right hand position. This trend is similar to 
the LTE models in Fig. 0]where systems with high ac- 
cretion rates are positioned in the upper left hand part 
of the Ha/H/3 versus H7/H/3 (abscissa versus ordinate) 
box and the lower mass flow rates are represented in 
the lower right hand location. The range of predicted 
line ratios are, however, wider for this model especially 
for Ha/H/3. In particular, the model predicts signifi- 
cantly non-flat Balmer decrements. Additionally, for 
disc models with temperatures between 15,000 K and 
8000 K, the predicted EW(/3) is much larger than for 
the optically thin LTE models. The predicted EW(/3) 
varies strongly as a function of No, going from ~ 200 
A for log/Vo = 11.5 to ~ 10 A for log/V = 14.0. 

3.2.4 Chromospheric model 

The prediction of Balmer line absolute strengths from 
disc chromospheres is made difficult because of the 
strong dependence of the un derlying continu um contri- 
butions on the disc model (|Fergusonl ll997h. However 
models of disc chromospheres have been made because 
of the probable existence of such chromospheres. Sit- 
uations where Q/k is a decreasing function of tem- 
perature, where Q is the energy generation and re is 
the Rosseland mean absorption coefficient in an atmo- 
sphere, can lead to an instability and a discontinuity in 
temperature in which t he temperature jumps h igher as 
the pressure decreases (|Shaviv fc Wehrsejl986t) . These 
are situations that lead to the formation of chromo- 
spheres and/or coronas. Such a situation was shown 
to exist in a pres cription accretion disc models by 
lAdam et all il988f) . Given th e prob able existence of 
disc chromospheres, IWilliamsl (|l99Rt l modelled emis- 
sion from chromospheres over a-disc models using ra- 
diative transfer methods, particularly computations us- 
ing the Feautrier method, similar to those used in the 
non-LTE isothermal models described in H3.2.3l above. 
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Figure 4: Photoi onization and LTE optically thin accretion disc models. The photoionization emission 
models are from lOsterbrockl l|1989ll and both Case A and Case B models (see text) are plotted. The 
Balmer emission line ratios from the LTE radiative transfer calculations assume an a-disc model with 
various mass transfer r ates ljWilliamsHl980|) . M is given in M Q yr 
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Figure 5: Scatterplots of Ha/H/3 versus H7/H/3 values from the non-LTE calculations 
Hydrogen number density, No, units are cm -3 ; i is disc inclination. 
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The B almer line emission ratios predicted by I Williams! 
lll995T) chromosphere models are summarized in Fig.|S| 
The associated predicted EW(/3) are small and vary 
from 1.02 A in model 6 to 32.27 A in model 11. 
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Figure 6: Chromospheric models of IWilliamsl 
(1995). For each model number, from 5 to 13, two 
outer disc radii are considered, 10 10 and 2 x 10 10 
cm with the model position moving from the upper 
left to lower right as the radius is increased. The 
assumed disc mass flow rate is 10 14 g s _1 for mod- 
els 5, 7 and 8, 10 15 for models 6, 9 and 12, 10 13 5 
for models 10 and 11, and 10 15 5 for model 13. The 
boundary layer luminosity, which photoionizes the 
chromospheric gas, is taken as compatible with the 
disc mass flow rate, except for models 8, 9 and 11 
where the boundary layer luminosity is that cor- 
responding to 10 times the disc mass flow rate. 
Finally the a disc viscosity parameter is assumed 
to be 1.0 for all models except for models 7, 12 and 
13 where a is assumed to be 0.5, 10.0 and 100.0 
respectively. 



4 Discussion 

lEchevarrial (|l988T l has compiled emission line data for 
CVs from several sources and produced a number of 
relevant scatterplots but did not subject the data to 
any formal multivariate statistical analysis. Echevarria 
reported a good fit (unquantified) of H Balmer emis- 
sion line rati o scatterplot with an LT E gas slab model 
computed bv lDrake fc Ultrichl l)l980T) in which position 
along roughly the principal component direction found 
in our present analysis depended only on disc temper- 
ature. Sources positioned further from the Drake & 
Ultrich model line in the H Balmer emission line ra- 
tio scatterplot space were found by Echevarria to have 
narrow H/3 equivalent widths. Here we have taken H 



Balmer emission line data from a single survey ijWilliamsl 
1983) and performed PCA and discriminant function 
analysis. 

The two multivariate analysis methods that we 
used literally give orthogonal results as we have indi- 
cated in jJ5]would be expected if the discriminant func- 
tion were found to be statistically significant. The first 
principal component direction is roughly in line with a 
least squares line fit through the 4D scatterplot data; 
the direction of the least squares line, V4, is 29.8° from 
the first PCA direction, Vi . The discriminant function 
directions V2 and V3 , on the other hand, form angles 
of 95.6° and 96.6°, respectively, with the first PCA di- 
rection. Further, the discriminant function directions 
are nearly orthogonal to the second and third princi- 
pal component directions and is most nearly aligned 
with the fourth principal component direction (within 
~ 15° — there is no statistical or geometrical reason 
for the discriminant function direction to match any 
of the principal component directions, its direction is 
strongly dependent on the direction of the difference 
between the mean vectors of the groups). The dif- 
ference between the first principal component direc- 
tion and the discriminant direction suggests that the 
physical processes responsible for the majority of the 
variance of the data (82% of the overall variance is 
in the first principal component direction) is a differ- 
ent one than what separates dwarf nova from other 
CVs. The first principal component direction reflects 
the fact that the Balmer decrement changes in a sys- 
tematic way from one system to the next. We also note 
that the "decrements" become "increments" for many 
systems (cf. Fig.Qwhere it may be seen that many ra- 
tios are greater than one for H7/H/3 to He/H/3). None 
of the models considered here (or any model that we 
are aware of) predict such "increments" (cf. Figs.^HUJ- 

iPattersonl I1984T) shows a good correlation for low- 
mass interacting binaries between the orbital period 
P or b and the mean mass transfer rate (M) from the 
secondary to the primary, as determined by a variety 
of methods, and presents an empirical relation 

(M) = 5.1±1 x icr 10 P 4 3 - 2±a2 M yr" 1 , (4) 

where P4 is P rt>/4 with P or b given in days. Although 
mass transfer rates for individual systems are known 
to vary widely from the rate giv en in equ ation Q, 
we can apply it to the data set of Williams] l|1983h to 
see if any average driven trends can id entified. Using 
Porb as given by iDownes et al.l i200ltl , M was com- 
puted using equation J1J and the correlation between 
M and the first principal component data values de- 
termined; a correlation of —0.436 (p = 1 x 10 -5 , two- 
tailed) was found. More directly, the correlation be- 
tween log Porb and the first principal component data 
values was found to be -0.664 (p = 5 x 10" 13 , two- 
tailed) , a result which may by appreciated in Fig. 
The implication is that the majority of the variance 
in the data along the first principal component direc- 
tion is due to variance in mass transfer rate. However, 
from inspection of Fig. [7| it appears that this conclu- 
sion may be strongly influenced by the T CrB systems 
which are known to have high mass flow rates ([Warned 
1995). Removal of the T CrB systems only reduces the 
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correlation between the first principal component and 
log P orb to -0.452 (p = 8x 1CT 6 , two-tailed) so the in- 
fluence of the T CrB systems did not overwhelmingly 
bias the correlation. 

The ordinate of the graph shown in Fig. Qroughly 
follows the direction (within 29.8°) of the line along 
which the data points lie in the 4D emission line ra- 
tio scatterplot space. The LTE model, as shown in 
Fig. 0] predicts a mass flow rate that increases as the 
first principal component coordinate increases. This 
variation is incompatible with M increasing with P or b , 
as given by equation (0J for example, and with an in- 
crease of P OI .b with a decreasing first principal compo- 
nent coordinate as shown in Fig. |7| A possible reason 
for the opposing correlations between (M) and the first 
principal component coordinate seen in the data and 
between M and the first principal component seen in 
the LTE models is that the LTE models may only be 
applicable to individual systems as an approximation 
when M varies because of disc instability and cannot 
be generalized to CV population behaviour where (M) 
varies because of orbital period differences. 

Discriminant function analyses on the basis of con- 
trasting dwarf novae with all other CVs or of contrast- 
ing dwarf novae with magnetic CVs both lead to sim- 
ilar separations. The mean discriminant function val- 
ues for dwarf novae are statistically significantly dif- 
ferent from other CVs and, from inspection of Figs. [5] 
and|^] the variance of H Balmer emission line discrim- 
inant function values is wider for dwarf novae than 



for other CVs 1 . Since the discriminant function direc- 
tions are nearly orthogonal to the first principal com- 
ponent direction, the physical processes that lead to 
the segregation of dwarf novae from other CVs in the 
discriminant function direction appear not to be re- 
lated directly to the differences in mass transfer rate if 
equation is taken as true. Supporting this conclu- 
sion, we additionally found no significant correlation 
between (M), as given by equation @), and discrimi- 
nant function values (p > 0.75, two-tailed) or between 
log P or b and discriminant function values (two-tailed 
p > 0.75 with the marginal exception of a correlation 
of —0.213, with two-tailed p = 0.041, between logP or b 
and the Class 2 discriminant function values — but 
any kind of correction for multiple comparisons makes 
this correlation non-significant). 

Visual comparison of Fig.^to Fig. [IJsuggests that 
the main source of H Balmer emission from T CrB 
systems is from nebula-like photoionization plasmas. 
Such photoionization may occur in the nebulousity sur- 
rounding the binary star, left behind from the recur- 
rent nova activity. The scatterplot values for the four 
T CrB systems represented in the dataset are, how- 
ever, more widely spaced than the values given by the 
theoretical calculations. This increase in spacing may 
be caused by the underlying absorption spectrum of 
the giant star which might tend to reduce Ha/H/3. 
An exception is RS Oph which shows line ratio val- 
ues consistent with a low temperature (5000 K) Case 
A (optically thin) situation. 

Possible sources of H Balmer emission in polars in- 
clude the stream between the secondary and the (mag- 
netic) threading region, gas falling down the magnetic 
field line, the heated surface of the stream close to 
the accretion zone or a reprocessing of high energy ac- 
cretion radiation on the surface of the secondary star. 
Since polars possess no disc and all other CVs likely do 
possess accretion discs of some sort, one would think 
that difference would discriminate polars from all other 
CVs in the 4D scatterplot space examined here. But 
that is not the case. Dwarf nova discs, however, may 
be distinguished from other CV discs in that the disc 
is unstable, and they undergo thermal cycling which 
leads to outbursts. The differences in H Balmer emis- 
sion line ratios between polars and non-dwarf nova 
CVs is smaller than the differences between dwarf no- 
vae and all other CVs, at least in the direction of the 
discriminant functions found here. There therefore ap- 
pear to be physical processes occurring in the accretion 
discs of dwarf novae that are different from those hap- 
pening with the other CV types and that there is a 
wide variation in the manner in which those unique 
dwarf nova processes occur. The obvious interpreta- 
tion is that, since the state of the accretion discs in 
dwarf novae are continuously changing, our data in- 
cludes discs at random phases between their last out- 
burst aadthenextone. 

lEchevarrial (I1988T) speculated that variation normal 

x It should be noted that interstellar reddening is a pos- 
sible source of variance for all groups, but this should not 
affect the difference in variances either between dwarf no- 
vae and other CVs or between dwarf novae and magnetic 
CVs. 
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Figure 8: Correlation of the equ ivalent width of 
H/3, EW(/3), as given hv lWilliamsl l)l 0.83ft . with the 
value of the discriminant function based on Class 
2. Shown is the least squares line fit through the 
data; a fit having a correlation of 0.422. No EW(/3) 
values for AM CVn stars were given by Williams. 



to our principal component direction was due to vari- 
ation in the equivalent width of the H/3 line, EW(/3). 
To test the hypothesis that variation along the Class 2 
discriminant function was related to equivalent width, 
we computed the correlation of the Class 2 discrimi- 
nant f unction values with EW(/3), as given by Williams 
(1983), and found a correlation of r = 0.422 (p = 
1 x 10 -4 , two-tailed). This correlation is illustrated 
in Fig. |H1 This correlation, in turn, implies a wide 
variance in EW(/3) of dwarf novae in the dataset, rel- 
ative to other CVs, from the wide variance of their 
discriminant function values (see Fig. |2J . The corre- 
lation of EW(/3) with the first principal component 
value is insignificant at r — —0.016 (p = 0.891) indi- 
cating that the cause of EW(/3) is closely related to 
the caus e of variation in the dis c rimin ant function di- 
rection. iPatterson fc Raymond! l|l985ft show that the 
ratio of X-ray flux (0.2 - 4.0 keV) to optical flux, 
F x /F v , is positively correlated with EW(/3) for CVs 
with accretion disks. If this were true for all CVs then 
the discriminant function direction would correspond 
to the relative X-ray intensity of the accretion region 
on the white dwarf. However, this interpretation is 
inconsistent with observed 0.001 < F x /F v < 10 for 
dwar f novae and 3 < F x /F v < 300 for polars ((Warned 
1995) when the values of the discriminant function as 
shown in Fig. [5] or [3] are considered. In other words, 
polars tend to have smaller EW(/3) than dwarf no- 
vae, not larger widths as would be implied by a larger 



F x /F v . The physical processes behind the production 
of Balmer emission in polars is not a continuous varia- 
tion of the processes that produce Balmer emission in 
dwarf novae. 

The Doppler width of the H/3 line, DW(/3) , is weakly 
correlated with both the Class 2 discriminant function 
value (r = 0.325, p — 0.006) and the first principal 
component value (r = 0.377, p = 0.001). Therefore 
factors influencing the DW(/3) value, such as broaden- 
ing due to disc rotation, have no coherent effect on the 
variation of Balmer emission line ratios in the popu- 
lation of CVs. But the variation in radiative transfer 
processes that give rise to variation of EW(/3) is co- 
herently related to the variation of emission line ratios 
in our sample of CVs, specifically in the discriminant 
function direction in which dwarf nova may be dis- 
tinguished from all other types of CVs. Inspection of 
Fig. [8] reveals that the variation of dwarf nova EW(/3) 
values in the discriminant function direction is much 
larger that the variation of other CV EW(/3) values, 
which tend to be smaller. This behaviour, of course, 
is simply a reflection of the behaviour of discriminant 
function values illustrated in Fig. [5] Since all dwarf 
nova contain unstable accretion discs ( Can nizzo et alJ 
1988) as a primary feature distinguishing them from 
other CVs, we may speculate that the variation of 
dwarf nova EW(/3) values is due to the variation of 
physical processes in the disc, perhaps related to ther- 
ma l cycling in the di sc as mentioned earlier. 

lEchevarrial l|1988ft essentially found variation along 
what we have defined as the first principal compo- 
nent direction to be a function of temperature while 
we hypothesized that it is essentially a function of 
mass transfer rate. The two conclusions are consis- 
tent because higher mass flow ra tes implies higher tem- 
peratures in the a-disc model flFrank. King fc Rainel 
(1995)). Echevarria found deviation in H Balmer line 
scatterplot space away from the LTE model of Drake 
& Ultrich to be related to the equivalent width of the 
H/3 line. We analogously find, more rigorously, that 
variation in the Class 2 discriminant function direc- 
tion, which is roughly orthogonal to the first principal 
component direction, is correlated with the equivalent 
width of the H/3 line. 



5 Conclusions 

Our analysis of hydrogen Balmer lines in CVs shows 
that the source of variation seen in the principal com- 
ponents analysis are correlated with the binary or- 
bital period, and the source of variation seen in the 
discriminant function analysis are correlated with the 
equivalent width of the H/3 line. Comparing models to 
the data produces one positive result. The ratios of 
hydrogen Balmer lines predicted by the photoionized 
model match the observed line ratios of the T CrB 
systems. In particular, the emission from RS Oph is 
consistent with photoionization from a ~5000 K op- 
tically thin (Case A) nebula. These systems, with a 
massive white dwarf and likely progenitors of type la 
supernovae, are hypothesized to be embedded i n a neb- 
ula produced by the frequent novae outbursts dWarneil 
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1995). Our analysis supports this scenario, as the ob- 
served Balmer emissions could originate from a pho- 
toionized nebula. Otherwise, the theoretical models 
considered in this work do not produce the wide range 
of line ratios along the first principal component direc- 
tion of the CV dataset. In particular, if we assume that 
the mass-transfer rate and the period are related by 
equation Q, the model of LTE emission from isother- 
mal medium predicts an opposite correlation between 
the observed first principal component direction and 
mass transfer rate than what is observed. 

Our analysis showed a wide variance of discrim- 
inant function values and a wide distribution of line 
ratios along the first principal component direction for 
dwarf novae. In light of the known thermal cycling 
in the accretion disc, regardless of whether the lines 
are emitted directly from the disc or are from repro- 
cessed radiation, both the Balmer line ratios and the 
line equivalent widths would be expected to vary with 
the dwarf nova outburst cycle. This may be a cause of 
the wide variance in the discriminant function values 
and of the wide distribution along the first principal 
component direction. As far as we are aware, system- 
atic observations of hydrogen Balmer emission lines 
throughout a dwarf nova cycle have not been reported. 

Our analysis here shows that the range of Balmer 
ratios seen along the first principle component direc- 
tion, including "inverted" Balmer decrements, can not 
be reproduced by the conventional models considered 
here. The more or less continuous variation of the 
line ratios from Balmer decrements to Balmer "incre- 
ments" for all CVs as a group prevents any clustering 
of the data that clearly points to a particular single 
underlying physical process. We therefore conclude 
that either the emission from different CV sub-classes 
have different locations (and hence line formation pro- 
cesses) within the system, or there are multiple regions 
for the hydrogen Balmer lines within a CV. The good 
matching between the data of T CrB systems and the 
photo-ionization nebular model supports the former. 
However, the latter is also likely, given that the UV 
radiation from the accreting white dwarf can ionize 
the cooler material in the outer accretion disc, the ac- 
cretion stream and the companion star. 
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